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Abstract 

We propose a novel stochastic method to exactly generate Brownian paths conditioned to start at 
an initial point and end at a given final point during a fixed time tf. These paths are weighted with a 
probability given by the overdamped Langevin dynamics. We show how these paths can be exactly 
generated by a local stochastic differential equation. The method is illustrated on the generation of 
Brownian bridges, Brownian meanders, Brownian excursions and constrained Ornstein-Uehlenbeck 
processes. In addition, we show how to solve this equation in the case of a general force acting on 
the particle. As an example, we show how to generate constrained path joining the two minima of a 
double-well. Our method allows to generate statistically independent paths, and is computationally 
very efficient. 
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I. INTRODUCTION 


Even after more than hundred years since its introduction, Brownian motion has remained 
a fundamental cornerstone of classical physics [T]. Simple Brownian motion and its variants 
have found numerous applications not just in natural sciences such as physics, chemistry, 
biology etc., but also in man made subjects such as finance, economics, queueing theory, 
search processes, and computer science amongst others (for reviews on the subject see [2:- 
E|). Some of these applications led to the studies of the statistical properties of simple 
variants of Brownian motion, generally referred to as ‘constrained Brownian motion’, i.e, 
a free Brownian motion conditioned to satisfy certain prescribed global constraints (see 
below for examples). More generally, one would also like to study other stochastic processes 
(thus going beyond the simple Brownian motion) in presence of one or more such global 
constraints. An interesting and challenging problem is to find simple and efficient algorithms 
that generate these constrained paths, for Brownian motion and other stochastic processes, 
with the correct statistical weight. Many of the concepts presented in this article through a 
physical presentation, can be found in the mathematical literature, for example in [7j. 

As a simple example, let us first start with a free Brownian motion B(t ) in one dimension, 
starting from the origin -B(O) = 0. The time evolution of B(t ) is governed by the simple 
Langevin equation 

dB . . 

dt-’M- (1) 

where 77 (f) is a Gaussian white noise with zero mean ( 77 (f)) = 0 and the correlator, 
(?7(t)77(t / )) = 2 DS(t — t') (D being the diffusion constant). It is easy to generate a trajec¬ 
tory of this process numerically, simply by simulating the time-discretized version of the 
Langevin equation 

B(t +At) = B{t) + V2DAff t , (2) 

where the noise £* = A/"(0,1) is a standard normally distributed random variable with zero 
mean and unit variance, drawn independently at each discrete time step. Three globally 
constrained variants (amongst others) of this simple 1-d Brownian motion that have been 
widely studied in the literature are the so called (i) Brownian bridge (ii) Brownian excursion 
and (iii) Brownian meander that we discuss next. 

Brownian bridge: A Brownian bridge xg(f) in 0 < f < tf is a Brownian motion, starting 
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Figure 1. (a) Two configurations of a Brownian bridge, starting at x B [t) = 0, and returning to 
x B (tf) = 0 at time tf (b) two configurations of a generalised Brownian bridge starting at x,b(0) = xq 
and ending at the final position XB(tf ) = Xf at time tf (xf fixed, but not necessarily zero). 

at xb( 0) = 0, that is constrained to come back to its starting point at a final time tf, 
i.e., x B {tf) = 0 (see Fig. 0 a)). A generalised Brownian bridge corresponds to those 
configurations of Brownian motions that start at #b(0) = Xq and end at x B (tf) — Xf (with 
Xf fixed, not necessarily the same as the starting position xq) at time tf (see Fig. 0b)). 

Brownian excursion: A Brownian excursion xg(i) in 0 < t < tf is a Brownian motion that 
starts at xe( 0) = 0, ends at x B (tf) = 0 (as in a bridge), but additionally is constrained 
to stay positive in between, i.e., for all 0 < t < tf (see Fig. (J2] a) for a typical excursion 
configuration). 

Brownian meander: A Brownian meander xm( 0 is a Brownian motion in [0, £/], that starts 
at Xm{ 0) = 0, is constrained to stay positive in the interval t £ [0 ,tf], but its final position 
x m( tf) is free, i.e., can be any positive number. For a typical meander configuration, see 
Fig. (@b). 

A natural question is: how does one generate a constrained path such as a bridge, excur¬ 
sion or a meander numerically with the correct statistical weight? Let us start with the ex¬ 
ample of a simple Brownian bridge x B (t) over [0, tf] with the constraint £#(0) = x B (tf ) = 0: 
how does one generate such a bridge? A naive algorithm would be to generate all possible 
configurations of a free Brownian motion over [0, tf] via the discretized Langevin equation 
0 starting from the origin and retain only those that do return to the origin at tf. Evi¬ 
dently, such a naive algorithm is too wasteful and would not lead to correct statistics for a 
finite number of samples. For the Brownian bridge, there is a however simple way out. One 
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Figure 2. (a) A typical configuration (schematic) of a Brownian excursion in the time interval [0, tf], 
that starts at xe( 0) = 0, ends at XEitf) = 0, and stays positive in the whole interval t E [0, tf] 
and (b) A typical configuration (schematic) of a Brownian meander in the time interval [0, tf], that 
starts at xm{ 0) = 0, stays positive in the time period [0, tf], but its final position XMitf) can be 
any positive number. 

first generates a free Brownian path B(t) via ([2]) starting from .8(0) = 0 and from this path, 
one constructs a new path Xs{t) via the following simple construction 

x B (t) = B(t) - y B(t f ) t E [0, t f ] . (3) 

By construction x B (t) in ([3]) satisfies the bridge condition x B (0) = ^s(t/) = 0. In addition, 
by computing the propagator for the constructed process xg(i), it is easy to show that 
indeed Xs{t ) has the correct statistical weight of a Brownian bridge. 

Similarly, one can generate an excursion Xsit), satisfying x^(0) = Xsitf) = 0 and xg(t) 
remaining positive in between, from a bridge configuration Xs{t) via the so called Vervaat 
construction [8| which proceeds as follows. First, we construct a bridge Xs{t) from a free 
Brownian path using ([3]). Let t m denote the time at which the bridge xs(t) achieves its 
minimum value say x m ; n (see Fig. (J3])). We divide the full time interval [0,i/] into two 
subintervals [0,t m ] and [t m ,tf]. Next we take the portion of the bridge path xsit) in [0, t m ] 
and slide it in the forward time direction by an interval tf [as shown in Fig. (|3|] and paste it 
to the other side of the rest of Xs(t) in the second subinterval [t m ,tf\. Due to the periodicity 
of the bridge, the newly slided and pasted section of the bridge will join smoothly with the 
part over [t m ,tf]. Next, we shift the origin to x m ; n , and the resulting process over the new 
time-interval [t m ,tf + t m ] is an excursion over [0, t/]- In other words, from the bridge Xsit), 
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Figure 3. Generation of an excursion from a bridge via Vervaat’s construction. One starts with a 
bridge configuration x B {t) with x_b(0) = x B (tf) = 0 (left side). The bridge achieves its minimum 
value x m j n at time t m (shown by the dashed vertical line). Divide the interval [0 ,tf] into two 
subintervals [0,i m ] (blue) and [t m ,tf] (red). Slide the left (blue) subinterval in the forward time 
direction by an amount tf (shown by the (green) arrow) and join it to the right of the subinterval 
Now shift the x axis so that x m i n denotes the new origin. This newly constructed cut 
and paste bridge configuration over t m , t m + t/] is indeed an excursion xg(t) over the time interval 

[0 ,t f \. 


we construct the excursion x E (t) via the construction 

x E (t) = x B (t m + t)~ x min t e [0 ,t f ] (4) 

where by periodicity x B {t + tf) = x B (t). By construction in Q, ^(0) = x B (tf) = 0 and 
> 0 for all 0 < t < tf. Moreover, Vervaat proved rigorously that the construction in 
Q indeed generates a Brownian excursion with the correct statistical weight [5J. 

Another simpler way to generate a Brownian excursion is to consider three independent 
Brownian bridges x B ^(t), x B ^{t) and x Bt3 (t) and then construct a new process 

x E (t ) = \J[x B ,i(t)] 2 + [x B< 2 (t)} 2 + [x B}3 (t)] 2 , (5) 

which just represents the radial part of a 3-dimensional Brownian bridge (3-d Bessel bridge). 
This new process can be rigorously shown to have the correct statistical weight of a one di¬ 
mensional Brownian excursion PSof. Similarly, a Brownian meander Xm{1) hr one dimen¬ 
sion, starting at %(0) = 0 and ending at Xuitf) = Xf can be generated from a 3-dimensional 
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Bessel bridge by the following construction PEO] 


x M (t) 


[xB,l(t)} 2 + [x B ,2(t)} 2 + [X B ,3 (t) + — XfY 


( 6 ) 


For a nice review on the connection between bridges, excursions and meanders, see mm 
and for recent applications of such constrained paths, see PESES]. 

These clever transformations connecting bridges, excursions and meanders to free Brow¬ 
nian motion are, of course, very efficient to generate them numerically. However, these 
transformations rely on the specific properties of a free Brownian motion and will, in gen¬ 
eral, not hold for a more general process such as a particle diffusing in an external potential. 
Moreover, one may want to generate paths which are not simply bridges, excursions or me¬ 
anders, but constrained in some other way. For example, in many applications in computer 
science, one needs to generate, say a Brownian excursion xsit) over t € [0 ,tf] with a fixed 
area under the excursion A = Jy f XE^t) dt p. fCTl fTTl [2QH22]. Hence, it would be much more 
desirable to build a recipe to construct an effective Langevin equation like (|2]) which (i) will 
automatically take into account the global constraints and (ii) will hold for more general 
stochastic processes such as a particle diffusing in an external potential. Once we have such 
a recipe for an effective Langevin equation, it can subsequently be easily time-discretized 
leading to an efficient method to generate constrained paths. The purpose of this paper is 
to provide precisely such a recipe and demonstrate it with several examples. Some of the re¬ 
sults to be discussed in this paper were probably known, at least partially, in the probability 
literature but in a language perhaps not easily accessible to physicists. One of the purposes 
of this paper is to unveil the recipe for an effective Langevin equation in a physicist friendly 
language. For a recent review on some aspects of constrained stochastic processes, namely, 
constraints on large deviations of a stochastic variable, see [23]. 

The rest of the paper is organized as follows: we first derive the effective Langevin 
equation which satisfies the boundary constraint. We then illustrate the method on four 
analytically soluble cases, namely the Brownian bridges, Brownian meanders, Brownian 
excursions and the Ornstein-LIhlenbeck process. We then show how this equation can be 
solved exactly numerically in the case of any landscape potential U{x ). 
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II. DERIVATION OF EFFECTIVE LANGEVIN EQUATION 

From now on, we assume that the system is driven by a force F(x, t ) and is subject to 
stochastic dynamics in the form of an overdamped Langevin equation. 

For the sake of simplicity, we illustrate the method on a one-dimensional system, the 
generalization to higher dimensions or larger number of degrees of freedom being straight¬ 
forward. We follow closely the presentation given in Ref. pa. 

The over damped Langevin equation reads 

dr 1 

- = -F(x(t),t)+ 77 (f) (7) 

where x(t) is the position of the particle at time f, driven by the force F(x, f), 7 is the friction 
coefficient, related to the diffusion constant D through the Einstein relation D = k B T/ 7 , 
where ks is the Boltzmann constant and T the temperature of the thermostat. In addition, 
77 (f) is a Gaussian white noise with moments given by 


< 77 (f) >= 0 

< >= ^ ~ f') 
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( 8 ) 

(9) 


The probability distribution P(x, t) for the particle to be at point x at time f satisfies a 
Fokker-Planck equation 




( 10 ) 


where /3 = 1 /k B T is the inverse temperature. This equation is to be supplemented by the 
initial condition P(ic,0) = S(x — x 0 ), where the particle is assumed to be at x 0 at time f = 0. 
To emphasize this initial condition, we will often use the notation P(x,t) = P(x, t\xo, 0). 

We now study the probability over all paths starting at Xo at time 0 and conditioned to 
end at a given point Xf at time tf, to fold the particle at point x at time f G [0, f/] - This is 
the generalized bridge. This probability can be written as 

1 


V{x,t) = 


where we use the notation 


P(x f ,t f |x 0 ,0) 
P(x,t ) = P(x, t\x 0 , 0) 


Q{x,t)P{x,t) 
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Q(x,t ) = P(xf,tf\x,t ) 


Indeed, the probability for a path starting from (x 0 , 0) and ending at ( Xf,tf ) to go through 
x at time t is the product of the probability P(x,t\x o ,0) to start at (x 0 ,0) and to end at 
(x,t) by the probability P(xf,tf\x,t) to start at (x, t) and to end at ( Xf,tf ). 

The equation satisfied by P is the Fokker-Planck equation mentioned above © , whereas 
that for Q is the so-called reverse or adjoint Fokker-Planck equation [26] given by 


d 4 =- D , 8 f& 

dt dx 2 dx 


( 11 ) 


It can be easily checked that the conditional probability V(x,t) satishes a new Fokker- 
Planck equation 


pp _ D &_ (&P _ 
dt dx V dx 


P F + 2 d -^)v 
dx J 


Comparing this equation with the initial Fokker-Planck (10) and Langevin ([7]) equations, 
one sees that it can be obtained from a Langevin equation with an additional potential force 


I- 


( 12 ) 


This equation has been originally obtained by Doob |27] in the probability literature and 
is known as the Doob transform of the Langevin equation (JT]) . It was also presented in 
the physics literature in the context of barrier crossing m It provides a simple recipe to 
construct a generalized bridge. It generates Brownian paths, starting at (ay, 0) conditioned 
to end at (xf,tf), with unbiased statistics. It is the additional term 2 D 91 ^ in the Langevin 
equation that guarantees that the trajectories starting at (a?o,0) and ending at (xf,tf) are 
statistically unbiased. Note that this additional term is a priori time dependent, and thus 
the effective force which conditions the paths is space and time dependent. 

In the following, we will specialize to the case where the force F is derived from a potential 
U(x). The bridge equation becomes 

dx D dU „ <9 In Q 


2D- 


+ v(t) 


(13) 


dt k B T dx dx 

In that case, the Fokker-Planck equation can be recast into an imaginary time Schrodinger 
equation [26], and the probability distribution function P can be written as 


Q(x,t) = P(x f ,t f \x,t ) = e -mv{x f )-u{x)) K Xf \ e -H(t f -t)\ x > 


(14) 









where the "quantum Hamiltonian" H is defined by 


H = —D 


d*_ 

dx 2 


DV(x) 


(15) 


and the potential V by 




2 dx J 2 dx 2 

We denote by M the matrix element of the Euclidian Schrodinger evolution operator 


(16) 


M(x,t)=<Xf\e H ( tf ^\x> 


Using eq.(14) for Q, one can write equation (13) as 


dx k B T d , . . , 

— = 2 -—In <xAe (/ ^ H \x>+n(t) 

dt 7 dx 


(17) 


We see on the above form that when t —> tf, the matrix element M{x,t) converges to 
5(xf — x), and it is this singular attractive potential which drives all the paths to Xf at time 

tf- 

Note that for large time tf, M{x,t ) is dominated by the ground state of H, namely 


<Xf\e t)JI \x>~e E ° ( ' tf f)v h 0 (x/)'ho(^) 


(18) 


It is well known that the ground state wave function of H is ^((e) = where 

J\f is the normalization constant. Indeed, it is easily checked that H |^ 0 ) = 0 which implies 
that d'o is an eigenstate of H with ground state energy E 0 = 0. Inserting these relations in 
eq. one recovers, as expected, the original unconditioned Langevin equation for large 


time tf. Note that although relation (18) holds only for a normalizable ground state wave 
function v I / o((e), it still yields the correct Langevin equation ([7]) for non-normalizable cases 
(such as the free Brownian motion U(x)= 0, for which Eq = 0 and H7 0 = 1/a/U where V is 
the volume). 

In order to build a representative sample of paths starting at (x 0 , 0) and ending at (. Xf,tf ), 


one must simply solve equation (13) for many different realizations of the random noise. Only 


the initial boundary condition is to be imposed, as the singular term in the equation imposes 
the correct final boundary condition. An important point to note is that all the trajectories 
generated by this procedure are statistically independent. 
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At this stage, we see that all the difficulty lies in calculating either the Green’s function 
Q(x,t) or the matrix element M(x,t). In general, this is not possible analytically except 
for a few special cases. In ref. |24] , a short time approximation to compute M(x,t) was 
presented, but in the following, we specialize on a few exactly analytically solvable cases. 


III. THE FREE BROWNIAN BRIDGE 


The simplest case is the case of a free particle in unrestricted space. In that case, the 
Green’s function Q(x,t ) can be easily computed and yields 


QOM) = 


(■ x f~ x r 
4 D(tf-t) 


sjATiD(tf — t) 

The corresponding conditioned Langevin equation becomes 

dx Xf — x 


dt tf — t 


+ v(t) 


(19) 


This equation (19) can be found in [7j j23]. It is completely universal as the drift term 
does not depend on the diffusion constant D. We show in Fig.Q a set of 500 trajectories 
starting at x$ = — 1 at time 0 and ending at Xf — +1 at time tf — 1, obtained for different 
noise histories. The time step used in the discretization is dt = 0.001. All these trajectories 
are statistically independent. The thick black curve is the mean trajectory. As can be seen 


from eq.(19), it is a straight line. 


IV. THE BROWNIAN MEANDER AND EXCURSION 

The Brownian meander is a Brownian path of a free particle in a half-space. In other 
words, the Brownian particle cannot penetrate in one of the half-planes. In Id, if we define 
the half-space as x > 0, the particle cannot go into the x < 0 subspace. Here the meander 
x M {t) starts, say at the origin = 0) = xq > 0, ends at XM{tf) — Xf > 0, while staying 
positive in t e [0,f/]. For the special case where both the initial and the final position tend 
to zero, the meander is called an excursion Xe(t). These constrained Brownian paths have 
been shown to be relevant to the anomalous diffusion of cold atoms in optical lattices |?5| . 

In the case of a general meander, the Green’s function Q(x, t) is that of a particle restricted 
to a half space. This Green’s function is well known and can be calculated by using the 
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Figure 4. A sample of 500 Brownian bridges 


method of images. In the case when the end point is fixed at Xf at time tf we obtain 
Q(x,t) = 


1 / _ 1 ( x f-x ) 2 _ 1 (xf+x ) 2 


4 D tf—t 


4 D tf—t 


y/4irD(t f -t) 

If the end point Xf at tf can be anywhere in the 1/2 plane Xf > 0, the Green’s function 
above has to be further integrated for Xf G [0, +oo] to obtain 


Q M (x,t) = / dxf Q(x,t) 

Jo 


= erf 


x 


( 20 ) 


where erf(x) = f* e _ “ 2 du. The corresponding Langevin equation for the meander then 


reads 


dx 


exp 


4 D(tf—t) 


+ !)(*) 


( 21 ) 


dt ys D(t, - t) erf( VjD ; f _„ ) 

The case of a Brownian excursion, where the extremity Xf is fixed, is generated by the 
Langevin equation 


dx 

dt 


1 (xf-x ) 2 

X f~ X ^ g _ 4D tf-t 
tf — t 


1 (xf +x ) 2 

+ ( TfB } tf ~ l 


1 (xf-x ) 2 

4 D tf—t 


l ( Xf+x ) 2 

4 D tf—t 


+ v(t) 


( 22 ) 
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Figure 5. A sample of 500 meanders 


Taking Xf —» 0, one obtains the effective Langevin equation for an excursion 


dx 2D f 
dt x \ 

This equation can be found in ref. [7]. 


x 


2 D(t f - t ) 


v(t) 


(23) 


Again, these two equations (21) and (23) can be solved easily by discretizing time. We 


show in Fig. (J5]) a set of 500 statistically independent meanders starting at Xq = 5 at time 0 
and ending anywhere in the upper half plane at time tf — 1. Similarly, we show in Fig.(j6]) 
a set of 500 statistically independent excursions starting at Xq = 0.01 at time 0 and ending 
at Xf = 0 at time tf — 1. In all cases, the independent trajectories are obtained for different 


noise histories and time step dt = 0.001. Note that generating an excursion by Eq. (23) is 


evidently much simpler than generating it by the Vervaat construction [S] discussed in the 
introduction. 

V. THE HARMONIC OSCILLATOR: THE ORNSTEIN-UHLENBECK PROCESS 


Another solvable case is that of the harmonic oscillator U(x) = Kx 2 /2. The uncondi¬ 
tioned Langevin equation is 

dnr 

(24) 


— = —Df3Kx + r)(t) 
dt 
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Figure 6. A sample of 500 excursions 

This is a so-called Ornstein-Uhlenbeck process |26j . The matrix element M(x. t) associated 
to this potential can be computed easily and one obtains 


M(x, t) = i exp f - 45illh ^ T) ((*? + J) cosh (DKfh) - 2xx f j 


(25) 


where t — tf — t and AT is a normalization constant. 

The associated conditioned Langevin equation is thus 




(26) 


which was obtained in [23] . 

Note that this equation is invariant when changing K into — K and thus we have this 
surprising result that paths going up a barrier are statistically identical to those going down 
a barrier. To illustrate this apparent paradox, we show in Fig. (Jr]) a set of 500 trajectories 
starting at Xo = — 1 at time 0 and ending at Xf — 0 at time tf — 1, obtained for different 
noise histories with dt = 0.001. The simulation was performed at temperature T = 0.1 with 
K = 1. 
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Figure 7. A sample of 500 Ornstein-Uhlenbeck trajectories starting at -1 and ending at 0 

VI. THE GENERAL CASE 


There are very few cases which are analytically solvable and we thus now discuss the 
generic non-solvable case. In order to compute the central quantity M(x,t ), we perforin a 
spectral expansion |28j 

M(x,t ) = 1 a;) 

= J2e- E ^~ t) 'i> a (xf)y a (x) (27) 

a 

where is the eigenstate of H with eigenvalue E a 


H^ a (x) = E a ^ { x) 


Note that although this spectral expansion is in principle valid only for normalizable 
eigenstates, it can easily be extended to the non-normalizable case by taking the limit of 
infinite volume. 


The bridge equation 0 becomes 


dx 

dt 


2D^~ _ f) 9x + n(t) 


(28) 


where we have used D = k B T/ 7 . 
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In order to be able to solve this equation numerically (by discretization for instance), 
one has to compute the eigenstates \h a and eigenvalues E a of the Hamiltonian H. In the 
case of low-dimensional system, this can be done very easily by diagonalizing the discretized 
Hamiltonian H which turns out to be a tridiagonal operator. 

VII. THE QUARTIC DOUBLE-WELL 

We illustrate the above method on the example of barrier crossing in Id (quartic poten¬ 
tial). 


U(x) = -(x‘ 2 - l) 2 

This potential has two minima at x — ±1, separated by a barrier of height 1/4. Note 
that V(x ) = {pV1 2) 2 — j3U"/2 is much steeper than U(x) and thus more confining around 
its minima. At low temperature, the potential V(x) has two minima at points close to ±1 
and one minimum at x — 0. 

The ground state of the Hamiltonian is To(x) ~ exp(— /3U(x)/2. The Hamiltonian is 
diagonalized by discretizing space and writing it in the form of a tridiagonal matrix. The 
diagonalization is performed using programs specific to tridiagonal matrices. All the exam¬ 
ples were performed at low temperature T = 0.05, where the barrier height is equal to 5 in 
units of ksT and the Kramers relaxation time, given by the inverse of the smallest non-zero 
eigenvalue of H , is equal to tk = 362.934 . 

On £g|8j we present a long trajectory (tf = 1000) obtained by solving the unconditioned 
Langevin eq.Q for a particle starting at Xo = — 1 at time 0 . 

The general pattern is that of the particle staying in one of the wells for a long time, 
then crossing very rapidly into the other well and back and forth. If one is interested more 
precisely in the transition region, our method allows to perform the simulation of the rare 
crossing events between the two wells and generate a large sample of statistically independent 
paths in this transition region. 

In hgj9j we show a set of 500 trajectories starting at Xq = — 1 at time 0 and ending at 


Xf = +1 at time tf = 10, obtained for different noise histories with dt = 0.001. In fig, 10 


we show a longer trajectory starting at Xq = —1 at time 0 and ending at Xf = +1 at time 
tf = 100 with dt = 0.001. 
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Figure 8. A sample long unconstrained Langevin trajectory in the double-well. 



Figure 9. A sample of 500 trajectories joining the two minima of the double-well. 


VIII. CONCLUSION 


We have presented in this paper a novel method in Physics to generate paths satisfying 
the Langevin overdamped dynamics, starting from an initial point and conditioned to end 
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Figure 10. A sample long trajectory in the double-well joining the two minima. 

at a given final point at a given time. One of the great advantages of this method is 
that trajectories are generated by a modified Langevin equation and are all statistically 
independent. In this paper, we emphasized applications to analytically solvable models, 
such as bridges, meanders, excursions, Ornstein-Uehlcnbeck processes, etc. In a forthcoming 
work, we will show how to generate paths with additional constraints, such as fixed area 
under the curve, fixed dissipation, etc. 
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